Soil moisture is an important boundary condition of earth-atmosphere exchanges, and it has been defined as an essential climate variable by GCOS. Vegetation optical depth is a physical variable to measure the attenuation of vegetation in microwave radiative transfer model, and it has been proved to be a good indicator of vegetation water content and biomass. This dataset uses the multi-channel collaborative algorithm (MCCA) to retrieve both soil moisture and polarized vegetation optical depth with SMAP brightness temperature. The algorithm uses a self-constraint relationship between land parameters and an analytical relationship between brightness temperature at different channels to perform the retrieval process. The MCCA does not depend on other auxiliary data on vegetation properties and can be applied to a variety of satellites. The soil moisture product from this dataset includes the soil moisture content in the unfrozen period and the liquid water content in the frozen period. Both horizontal- and vertical-polarization vegetation optical depth are retrieved. So far as we know, it is the first polarization-dependent vegetation optical depth product at L-band. This dataset was validated by 19 dense soil moisture observation networks (9 core validation sites used by SMAP team and 13 sites not used by them), and the widely used soil climate analysis network (SCAN). It was found that ubRMSE (unbiased root mean square error) of MCCA retrieved soil moisture is generally smaller than that of other SMAP products.
ZHAO Tianjie, PENG Zhiqing , YAO Panpan, SHI Jiancheng
Based on the data of GF-1 and GF-2 in China, the freeze-thaw disaster distribution data of Qinghai Tibet project corridor is produced by using the deep learning classification method and manual visual interpretation and correction. The geographical range of the data is 40km along the Xidatan Anduo section of Qinghai Tibet highway. The data include the distribution data of thermokast lakes and the distribution data of thermal melting landslides. The dataset can provide data basis for the research of freeze-thaw disaster and engineering disaster prevention and reduction in Qinghai Tibet engineering corridor. The spatial distribution of freezing and thawing disasters within 40km along the Xidatan-Anduo section of Qinghai Tibet highway is self-made based on the domestic GF-2 image data. Firstly, the deep learning method is used to extract the mud flow terrace block from GF-2 data; Then, ArcGIS is used for manual editing.
NIU Fujun, LUO Jing LUO Jing
This phenological data is based on the MOD13A2 data of the Qinghai Tibet Plateau from 2000 to 2015 (with a temporal resolution of 16 days and a spatial resolution of 1km). The NDVI curve is fitted using the segmented Gaussian function in the TIMESAT software. The spring phenology, autumn phenology and the length of the growth season are extracted using the dynamic threshold method. The thresholds of spring phenology and autumn phenology are set to 0.2 and 0.7 respectively. The phenological data were masked. Among them, the mask rules are: 1) The maximum value of NDVI must be met between June and September; 2) The average value of NDVI from June to September shall not be less than 0.2; 3) The average NDVI in winter shall not exceed 0.3.
ZU Jiaxing , ZHANG Yangjian
The normalized difference vegetation index (NDVI) can accurately reflect the surface vegetation coverage. At present, NDVI time series data based on spot / vegetation and MODIS satellite remote sensing images have been widely used in the research of vegetation dynamic change monitoring, land use / cover change detection, macro vegetation cover classification and net primary productivity estimation at various scales. Evi is similar to the normalized difference vegetation index (NDVI) and can be used to quantify vegetation greenness. However, evi corrects for some atmospheric conditions and canopy background noise and is more sensitive in areas with dense vegetation. It contains an "L" value to adjust the canopy background, a "C" value as the atmospheric drag coefficient, and a value from the blue band (b). These enhancements allow the ratio between R and NIR values to be calculated exponentially while reducing background noise, atmospheric noise and saturation in most cases. This research work mainly focuses on post-processing NDVI and evi data, and gives a more reliable vegetation situation of the Qinghai Tibet Plateau in 2013 and 2018 through transformation of projection coordinate system, data fusion, maximum value synthesis method, elimination of outliers and clipping. The spatial resolution of the data is 0.05 °, and the temporal resolution is month.
YE Aizhong
NDVI reflects the background effects of plant canopy, such as soil, wet ground, snow, dead leaves, roughness, etc., and is related to vegetation cover. It is one of the important parameters to reflect the crop growth and nutrient information. According to this parameter, the N demand of crops in different seasons can be known, which is an important guide to the reasonable application of N fertilizer. Correct NDVI (C-NDVI) is the value of NDVI after excluding the influence of climate elements (temperature, precipitation, etc.) on NDVI. Taking precipitation as an example, studies on the lag effect of precipitation on vegetation growth show that the lag time of precipitation effects varies in different regions due to differences in vegetation composition and soil types. In this study, we post-processed the MODIS NDVI data and firstly correlated the NDVI value of the current month with the precipitation of the current month, the average value of the precipitation of the current month with that of the previous month, and the average value of the precipitation of the current month with that of the previous two months to determine the optimal lag time. The NDVI was regressed on precipitation and air temperature to obtain the correlation coefficients, and then the corrected NDVI values were calculated by the difference between the MODIS NDVI and the NDVI regressed on climate factors. We corrected NDVI using climate data to give reliable vegetation correction indices for the circum-Arctic Circle (range north of 66°N) and the Tibetan Plateau (range 26°N to 39.85°N and 73.45°E to 104.65°E) for 2013 and 2018. The spatial resolution of the data is 0.5 degrees and the temporal resolution is monthly values.
YE Aizhong
Data content: this data set is the historical archived satellite data of the domestic high score series (GF1 / 2 / 3 / 4) in the key river and lake research areas of the Qinghai Tibet Plateau from 2015 to 2020, which can cover the typical river and lake areas for effective monitoring. The time range of the data is from 2015 to 2020. Data source and processing method: the data are level 1 products. After equalizing radiation correction, the changes affecting the sensors are corrected by the equalizing functions of different detectors. Some data are based on the Landsat 8 images in the same period as the base map, and control points are selected for geometric correction of the images. Then, orthophoto correction is carried out based on DEM data, and band fusion processing is carried out for the corresponding data. Data quality description: the Gaofen series satellites are processed by the China Resources Satellite Application Center. There are raw data received by the satellite ground receiving station of the Chinese Academy of Sciences and processed products at all levels. Among them, level 1a (pre-processing level radiometric correction image product): image data processed by data analysis, uniform radiometric correction, noise removal, MTFC, CCD splicing, band registration, etc; And provide RPC files for satellite direct attitude orbit data production. Refer to the data website of China Resources Satellite Application Center for details. Data application achievements and prospects: the data are domestic high-resolution data with high resolution, which can be used to monitor the changes of the Qinghai Tibet Plateau as a water tower in Asia and the generated images, and test the accuracy of other data in the region
QIU Yubao
The dataset is the remote sensing image data ofGF-1 satellite in the Qinghai-Tibet engineering corridor obtained by China High Resolution Earth Observation Center. After the fusion processing of multispectral and panchromatic bands, the image data with a spatial resolution of 2 m is obtained. In the process of obtaining ground vegetation information, the classification technology of combining object-oriented computer automatic interpretation and manual interpretation is adopted, The object-oriented classification technology is to collect adjacent pixels as objects to identify the spectral elements of interest, make full use of high-resolution panchromatic and multispectral data space, texture and spectral information to segment and classify, and output high-precision classification results or vectors. In actual operation, the image is automatically extracted by eCognition software. The main processes are image segmentation, information extraction and accuracy evaluation. After verification with the field survey, the overall extraction accuracy is more than 90%.
NIU Fujun
Fractional Vegetation Cover (FVC) refers to the percentage of the vertical projected area of vegetation to the total area of the study area. It is an important indicator to measure the effectiveness of ecological protection and ecological restoration. It is widely used in the fields of climate, ecology, soil erosion and so on. FVC is not only an ideal parameter to reflect the productivity of vegetation, but also can play a good role in evaluating topographic differences, climate change and regional ecological environment quality. This research work is mainly to post process two sets of glass FVC data, and give a more reliable vegetation coverage of the circumpolar Arctic Circle (north of 66 ° n) and the Qinghai Tibet Plateau (north of 26 ° n to 39.85 °, east longitude 73.45 ° to 104.65 °) in 2013 and 2018 through data fusion, elimination of outliers and clipping.
YE Aizhong
Photosynthetically active radiation (PAR) is fundamental physiological variable driving the process of material and energy exchange, and is indispensable for researches in ecological and agricultural fields. In this study, we produced a 35-year (1984-2018) high-resolution (3 h, 10 km) global grided PAR dataset with an effective physical-based PAR model. The main inputs were cloud optical depth from the latest International Satellite Cloud Climatology Project (ISCCP) H-series cloud products, the routine variables (water vapor, surface pressure and ozone) from the ERA5 reanalysis data, aerosol from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) products and albedo from Moderate Resolution Imaging Spectroradiometer (MODIS) product after 2000 and CLARRA-2 product before 2000. The grided PAR products were evaluated against surface observations measured at seven experimental stations of the SURFace RADiation budget network (SURFRAD), 42 experimental stations of the National Ecological Observatory Network (NEON), and 38 experimental stations of the Chinese Ecosystem Research Network (CERN). The instantaneous PAR was validated at the SURFRAD and NEON, and the mean bias errors (MBEs) and root mean square errors (RMSEs) are 5.6 W m-2 and 44.3 W m-2, and 5.9 W m-2 and 45.5 W m-2, respectively, and correlation coefficients (R) are both 0.94 at 10 km scale. When averaged to 30 km, the errors were obviously reduced with RMSEs decreasing to 36.3 W m-2 and 36.3 W m-2 and R both increasing to 0.96. The daily PAR was validated at the SURFRAD, NEON and CERN, and the RMSEs were 13.2 W m-2, 13.1 W m-2 and 19.6 W m-2, respectively at 10 km scale. The RMSEs were slightly reduced to 11.2 W m-2, 11.6 W m-2, and 18.6 W m-2 when upscaled to 30 km. Comparison with the other well-known global satellite-based PAR product of the Earth's Radiant Energy System (CERES) reveals that our PAR product was a more accurate dataset with higher resolution than the CRERS. Our grided PAR dataset would contribute to the ecological simulation and food yield assessment in the future.
TANG Wenjun
Global land cover data are key sources of information for understanding the complex interactions between human activities and global change. FROM-GLC (Finer Resolution Observation and Monitoring of Global Land Cover) from Tsinghua is the 30 m resolution global land cover maps produced. The Global land cover data of all 34 key nodes of pan-third pole region are produced through analyse by argis. The classfication system is crop(10), forest(20), grass(30), shrbu(40), wetland(50), water(60), tundra(70), impervious(80), Bareland(90), snow/ice(100), cloud(120). Finally, This data set serves as the research basis for all remote sensing data and provides baseline data for the project.
GE Yong, LING Feng, ZHANG Yihang
This data set is daily surface albedo product over Tibet plateau region from 2002 to 2020 with a spatial resolution of 0.00425°. The MODIS reflectance data product was used to retrieve the Extended Multi-Sensor Combined BRDF Inversion (EMCBI) Model which has coupled with topographic effects with assistance of a BRDF priori-knowledge. The daily BRDF was retrieved in a 5-day period to collect multi-angular information from MODIS observations. And then the daily albedo is estimated, where the black sky albedo was calculated at local noon. MODIS surface reflectance data (MOD09GA and MYD09GA) are downloaded from the official website. The albedo product is quality-controlled with better temporal and spatial continuity in Tibet plateau area. The validation results show that it meets the accuracy requirements of albedo application with higher precisions comparing to the other similar products. And thus, this product is useful for the long-term environmental monitoring and radiation energy budget research study.
YOU Dongqin, YOU Dongqin, TANG Yong, TANG Yong, TANG Yong, HAN Yuan HAN Yuan
The basic data set of remote sensing for ecological assets assessment of the Qinghai-Tibet Plateau includes the annual Fraction Vegetation Coverage (FVC), Net Primary Productivity (NPP) and Leaf Area Index (LAI) of the Qinghai-Tibet Plateau since 2000, and other ecological parameters based on remote sensing inversion. The FVC data are mainly developed from MODIS NDVI data. Based on pixel dichotomy model, the vegetation coverage model is developed by using multi-scale remote sensing images, combining with high precision remote sensing parameters such as vegetation community type and distribution characteristics, and the mixed pixel decomposition method is used to construct the vegetation coverage model. All data could be used only after the permission of the data distributor.
LIU Wenjun
This dataset is a pixel-based maximum fractional vegetation cover map within the Yellow River source region on the Qinghai-Tibet Plateau, with an area of about 44,000 square kilometers. Based on the time series images acquired from MODIS with a resolution of 250 m and Landsat-8 with a resolution of 30 m in 2015 during the vegetation growing season, the data are derived using dimidiate pixel model and time interpolation. The spatial resolution of the image is 30 m, using the WGS 1984 UTM projected coordinate system, and the data is in the format of grid.
WANG Guangjun
Based on 11 well-acknowledged global-scale microwave remote sensing-based surface soil moisture products, and with 9 main quality impact factors of microwave-based soil moisture retrieval incorporated, we developed the Remote Sensing-based global Surface Soil Moisture dataset (RSSSM, 2003~2020) through a complicated neural network approach. The spatial resolution of RSSSM is 0.1°, while the temporal resolution is approximately 10 days. The original dataset covered 2003~2018, but now it has been updated to 2020. RSSSM dataset is outstanding in terms of temporal continuity, and has full spatial coverage except for snow, ice and water bodies. The comparison against the global-scale in-situ soil moisture measurements indicates that RSSSM has a higher spatial and temporal accuracy than most of the frequently-used global/regional long-term surface soil moisture datasets. In addition, although RSSSM is remote sensing based, without the incorporation of any precipitation data or records, its interannual variation generally conforms with that of precipitation (e.g., the GPM IMERG precipitation data) and Standardized Precipitation Evapotranspiration Index (SPEI). Moreover, RSSSM can also reflect the impact of human activities, e.g., urbanization, cropland irrigation and afforestation on soil moisture changes to some degree. The data is in ‘Tiff’ format, and the size after compression is 2.48 GB. The relevant data describing paper has been published in the Journal ‘Earth System Science Data’ in 2021.
CHEN Yongzhe, FENG Xiaoming, FU Bojie
The global monthly all-sky land surface temperature (2000-2020) is produced by the method from Chen et al. 2017 JHM.
CHEN Xuelong, BOB Su, MA Yaoming
The data was obtained from the 30-second global elevation dataset developed by the US Geological Survey (USGS) and completed in 1996. Downloaded the data from the NCAR and UCAR Joint Data Download Center (https://rda.ucar.edu/datasets/ds758.0/) and redistributed it through this data center. GTOPO30 divides the world into 33 blocks. The sampling interval is 30 arc seconds, which is 0.00833333333333333 degrees. The coordinate reference is WGS84. The DEM is the distance from the sea level in the vertical direction, ie the altitude, in m, the altitude range from -407 to 8752, the ocean depth information is not included here, the negative value is the altitude of the continental shelf; the ocean is marked as -9999, the elevation above the coastline is at least 1; the island less than 1 square kilometer is not considered. In order to facilitate the user's convenience, on the basis of the block data, splice 10 blocks in -10S-90N and 20W-180E without any resampling processing. This data file is DEM_ptpe_Gtopo30.nc
HE Yongli
Thematic data on desertification (land desertification, salinization and vegetation degradation) in Central Asia, includes three parts: Distribution Map of Sandy Land in Central Asia, Distribution Map of Salinized Land in Central Asia and Distribution Map of Land Vegetation Degradation in Central Asia. The spatial resolution of the data is 1km, the time resolution is in 2015. The data produced by the key laboratory of remote sensing and GIS, Xinjiang institute of ecology and geography, Chinese Academy of Sciences. Data production Supported by the Strategic Priority Research Program of Chinese Academy of Sciences, Grant No. XDA20030101.
XU Wenqiang
Based on a large number of measured aboveground biomass data of grassland, the temperate grassland types were divided according to the vegetation type map of China in 1980s Based on the Landsat remote sensing data of engine platform, the random forest model of grassland aboveground biomass and remote sensing data was constructed for different grassland types. On the basis of reliable verification, the annual aboveground biomass of grassland from 1993 to 2019 was estimated, and the annual spatial data set of aboveground biomass of temperate grassland in Northern China from 1993 to 2019 was formed. Aboveground biomass is defined as the total amount of organic matter of vegetation living above the ground in unit area. The original grid value has been multiplied by a factor of 100, unit: 0.01 g / m2 (g / m2). This data set can provide a scientific basis for the dynamic monitoring and evaluation of temperate grassland resources and ecological environment in northern China.
ZHANG Na
The distribution data of Central Asia desert oil and gas fields are in the form of vector data in ". SHP". Including the distribution of oil and gas fields and major urban settlements in the five Central Asian countries. The data is extracted and cut from modis-mcd12q product. The spatial resolution of the product is 500 m, and the time resolution is 1 year. IGBP global vegetation classification scheme is adopted as the classification standard. The scheme is divided into 17 land cover types, among which the urban data uses the construction and urban land in the scheme. The data can provide data support for the assessment and prevention of sandstorm disasters in Central Asia desert oil and gas fields and green town.
GAO Xin
Thematic data on desertification in Western Asia, includes two parts: Distribution Map of Sandy Land in Western Asia, Distribution Map of Grassland Degradation in Western Asia. The spatial resolution of the data is 30m. The data produced by the key laboratory of remote sensing and GIS, Xinjiang institute of ecology and geography, Chinese Academy of Sciences, the spatial resolution of data is 30 m. Data production Supported by the Strategic Priority Research Program of Chinese Academy of Sciences, Grant No. XDA20030101. The map of artificial oasis pattern in Amu river basin is based on Landsat TM and ETM image data in 2015. Firstly, with the help of eCognition software, the object-oriented classification is carried out. Secondly, the classification results are checked and corrected manually.
Evapotranspiration over the Qinghai Tibet Plateau is calculated by etwatch, a land surface evapotranspiration remote sensing model based on multi-scale and multi-source data. Etwatch adopts the method of combining the residual term method with P-M formula to calculate evapotranspiration. Firstly, according to the characteristics of the data image, the suitable model is selected to retrieve the evapotranspiration on a sunny day; the remote sensing model is often lack of data because the weather conditions can not obtain a clear image. In order to obtain the daily continuous evapotranspiration, the penman Monteith formula is introduced, and the evapotranspiration results on a sunny day are regarded as the "key frame", and the surface impedance information of the key frame is used as the basis to construct the surface impedance Based on the daily meteorological data, the time series data of evapotranspiration are reconstructed. Through the data fusion model, the high spatial and temporal resolution evapotranspiration data set is constructed by combining the low and medium resolution evapotranspiration temporal variation information with the high resolution evapotranspiration spatial difference information, so as to generate the 8 km resolution evapotranspiration of the Qinghai Tibet Plateau Data sets (1990-2015).
WANG Xiaofeng
Sentine-1 SAR data were used to monitor the permafrost of Biuniugou in Heihe River Basin of Qinghai-Tibet Plateau. Based on the Sentine-1 SAR image of Bison Valley from 2014 to 2018, the active layer thickness in the study area was estimated by using the small baseline set time series InSAR (DSs-SBAS) frozen soil deformation monitoring method based on distributed radar target, combined with SAR backscattering coefficient, MODIS surface temperature and Stefan model. The results show that the thickness of active layer is between 0.8 m and 6.6 m, with an average of about 3.3 M. It is of great significance to carry out large-scale and high-resolution monitoring.
JIANG Liming
Data Set of Key Elements of Desertification in Typical Watershed of Central and Western Asia includes four parts: distribution and change of agricultural land of Amu River Basin, distribution and change of grassland of Amu River Basin, distribution and change of shrub land of Amu River Basin, distribution and change of forests of Amu River Basin. the spatial resolution of data is 30 m. All the data is based on Landsat TM/ETM image data in 1990, 2000 and 2010. The data produced by the key laboratory of remote sensing and GIS, Xinjiang institute of ecology and geography, Chinese Academy of Sciences. Data production Supported by the Strategic Priority Research Program of Chinese Academy of Sciences, Grant No. XDA20030101.
This dataset includes Fraction Vegetation Coverage (FVC) data for five countries in Central Asia (Kazakhstan, Kyrgyzstan, Tajikistan, Turkmenistan, and Uzbekistan) during 2010, 2015 and 2020. The data is calculated from the MODIS-NDVI data set (product number MOD13A2.006) based on the empirical relationship between FVC in arid areas and NDVI. The product has a time resolution of 1 year and a spatial resolution of 1 km. The algorithm selects the best available pixel value based on low cloud, low detection angle and highest NDVI value from all the observation data of the year, and performs conversion.
XU Xiaofan, TAN Minghong
Surface solar irradiance (SSI) is one of the products of FY-4A L2 quantitative inversion. It covers a full disk without projection, with a spatial resolution of 4km and a temporal resolution of 15min (there are 40 observation times in the whole day since 20180921, except for the observation of each hour, there is one observation every 3hr before and after the hour）, and the spectral range is 0.2µ m～5.0 µ m. The output elements of the product include total irradiance, direct irradiance on horizontal plane and scattered irradiance, the effective measurement ranges between 0-1500 w / m2. The qualitative improvement of FY-4A SSI products in coverage, spatial resolution, time continuity, output elements and other aspects makes it possible to further carry out its fine application in solar energy, agriculture, ecology, transportation and other professional meteorological services. The current research results show that the overall correlation of FY-4A SSI product in China is more than 0.75 compared with ground-based observation, which can be used for solar energy resource assessment in China.
SHEN Yanbo, HU Yueming, HU Xiuqing
This daily land surface kernel-driven BRDF model's coeciffients proudct is with a spatl resolution of 0.02 ° x 0.02 ° over the Tibet Plateau in 2016. Multi-sensory data is used to retrieve the the kernel-driven BRDF model and coupled with topographic effects, and prior knowledge is introduced for quality control inversion. The high-precision BRDF of good spatial-temporal continiuty is retrieved by combining MODIS reflectance data (a polar orbiting satellite) and himawari-8 AHI land surface reflectance (a geostationary satellite ). MODIS lans surface reflectance data and AHI TOA reflectance data are downloaded from the official websites. After registration, atmospheric correction and other processing, the daily resolution BRDF is synthesized with a period of 5 days. Compared with similar products, it has more advantages in capturing rapidly changing surface features, and has better temporal and spatial continuity with the shortest composition period. It can effectively support angular effects correction and the BRDF-releated parameters' retrieval.
WEN Jianguang , TANG Yong, TANG Yong, TANG Yong, YOU Dongqin YOU Dongqin
Based on the medium resolution long time series remote sensing image Landsat, the data set obtained six periods of ecosystem type distribution maps of the Qinghai Tibet Plateau in 1990 / 1995 / 2002 / 2005 / 2010 / 2015 through image fusion, remote sensing interpretation and data inversion, and made the original ecological base map of the Qinghai Tibet Plateau in 25 years (1990-2015). According to the area statistics of various ecosystems in the Qinghai Tibet Plateau, the area of woodland and grassland decreased slightly, the area of urban land, rural residential areas and other construction land increased, the area of rivers, lakes and other water bodies increased, and the area of permanent glacier snow decreased from 1990 to 2015. The atlas can be used for the planning, design and management of ecological projects in the Qinghai Tibet Plateau, and can be used as a benchmark for the current situation of the ecosystem, to clarify the temporal and spatial pattern of major ecological projects in the Qinghai Tibet Plateau, and to reveal the change rules and regional differences of the pattern and function of the ecosystem in the Qinghai Tibet Plateau.
ZHAO Hui, WANG Xiaodan
The data set mainly includes the ice observation frequency (ICO) of north temperate lakes in four periods from 1985 to 2020, as well as the location, area and elevation of the lakes. Among them, the four time periods are 1985-1998 (P1), 1999-2006 (P2), 2007-2014 (P3) and 2015-2020 (P4) respectively, in order to improve the "valid observation" times in the calculation period and improve the accuracy. The ICO of the four periods is calculated by the ratio of "icing" times and "valid observation" times counted by all Landsat images in each period. Other lake information corresponds to the HydroLAKEs data set through the "hylak_id" column in the table. In addition, the data only retains about 30000 lakes with an area of more than 1 square kilometer, which are valid for P1-P4 observation. The data set can reflect the response of Lake icing to climate change in recent decades.
WANG Xinchi
This dataset contains the glacier outlines in Qilian Mountain Area in 2020. The dataset was produced based on classical band ratio criterion and manual editing. Chinese GF series images collected in 2020 were used as basic data for glacier extraction. Google images and Map World images were employed as reference data for manual adjusting. The dataset was stored in SHP format and attached with the attributions of coordinates, glacier ID and glacier area. Consisting of 1 season, the dataset has a spatial resolution of 2 meters. The accuracy is about 1 pixel (±2 meter). The dataset directly reflects the glacier distribution within the Qilian Mountain in 2020, and can be used for quantitative estimation of glacier mass balance and the quantitative assessment of glacier change’s impact on basin runoff.
Li Jia Li Jia LI Jia LI Jia
These data contain two data files: GLOBELAND30 TILES (raw data) and TIBET_ GLOBELAND30_MOSAIC (mosaic data). The raw data were downloaded from the Global Land Cover Data website (GlobalLand3) (http://www.globallandcover.com) and cover the Tibetan Plateau and surrounding areas. The raw data were stored in frames, and for the convenience of using the data, we use Erdas software to splice and mosaic the raw data. The Global Land Cover Data (GlobalLand30) is the result of the “Global Land Cover Remote Sensing Mapping and Key Technology Research”, which is a key project of the National 863 Program. Using the American Landsat images (TM5, ETM+) and Chinese Environmental Disaster Reduction Satellite images (HJ-1), the data were extracted by a comprehensive method based on pixel classification-object extraction-knowledge checks. The data include 10 primary land cover types—cultivated land, forest, grassland, shrub, wetland, water body, tundra, man-made cover, bare land, glacier and permanent snow—without extracting secondary types. In terms of accuracy assessment, nine types and more than 150,000 test samples were evaluated. The overall accuracy of the GlobeLand30-2010 data is 80.33%. The Kappa indicator is 0.75. The GlobeLand30 data use the WGS84 coordinate system, UTM projection, and 6-degree banding, and the reference ellipsoid is the WGS 84 ellipsoid. According to different latitudes, the data are organized into two types of framing. In the regions of 60° north and south latitudes, the framing is carried out according to a size of 5° (latitude) × 6° (longitude); in the regions of 60° to 80° north and south latitudes, the framing is carried out according to a size of 5° (latitude) × 12° (longitude). The framing is projected according to the central meridian of the odd 6° band. GLOBELAND30 TILES: The original, unprocessed raw data are retained. TIBET_ GLOBELAND30_MOSAIC: The Erdas software is used to mosaic the raw data. The parameter settings use the default value of the raw data to retain the original, and the accuracy is consistent with that of the downloading site.
CHEN Jun
This data set includes the monthly average actual evapotranspiration of the Tibet Plateau from 2001 to 2018. The data set is based on the satellite remote sensing data (MODIS) and reanalysis meteorological data (CMFD), and is calculated by the surface energy balance system model (SEBS). In the process of calculating the turbulent flux, the sub-grid scale topography drag parameterization scheme is introduced to improve the simulation of sensible and latent heat fluxes. In addition, the evapotranspiration of the model is verified by the observation data of six turbulence flux stations on the Tibetan Plateau, which shows high accuracy. The data set can be used to study the characteristics of land-atmosphere interaction and the water cycle in the Tibetan Plateau.
HAN Cunbo, MA Yaoming, WANG Binbin, ZHONG Lei, MA Weiqiang*, CHEN Xuelong, SU Zhongbo
Nighttime light remote sensing has been an increasingly important proxy for human activities including socioeconomics and energy consumption. Defense Meteorological Satellite Program-Operational Linescan System from 1992 to 2013 and Suomi National Polar-Orbiting Partnership-Visible Infrared Imaging Radiometer Suite since 2012 are the most widely used datasets. Despite urgent needs for long-term products and pilot explorations in synthesizing them, the publicly available long-term products are limited. We propose a Night-Time Light convolutional Long Short-Term Memory (NTLSTM) network, and apply the network to produce annual Prolonged Artificial Nighttime-light DAtaset (PANDA) in China from 1984 to 2020. Model assessments between modelled and original images show that on average the Root Mean Squared-Error (RMSE) reaches 0.73, the coefficient of determination (R2) reaches 0.95, and the linear slope is 0.99 at pixel level, indicating a high confidential level of the data quality of the generated product. In urban areas, the modelled results can well capture temporal trends in newly built-up areas but slightly underestimate the intensity within old urban cores. Socioeconomic indicators (built-up areas, Gross Domestic Product, population) correlates better with the PANDA than with previous products in the literature, indicating its better potential in finding different controls of nighttime-light variances in different phases. Besides, the PANDA delineates different urban expansion types, outperforms other products in representing road networks, and provides potential nighttime-light sceneries in early years. PANDA provides the opportunity to better bridge the cooperation between human activity observations and socioeconomic or environmental fields
ZHANG Lixian, REN Zhehao, CHEN Bin, GONG Peng, FU Haohuan, XU Bing
The long-term evolution of lakes on the Tibetan Plateau (TP) could be observed from Landsat series of satellite data since the 1970s. However, the seasonal cycles of lakes on the TP have received little attention due to high cloud contamination of the commonly-used optical images. In this study, for the first time, the seasonal cycle of lakes on the TP were detected using Sentinel-1 Synthetic Aperture Radar (SAR) data with a high repeat cycle. A total of approximately 6000 Level-1 scenes were obtained that covered all large lakes (> 50 km2) in the study area. The images were extracted from stripmap (SM) and interferometric wide swath (IW) modes that had a pixel spacing of 40 m in the range and azimuth directions. The lake boundaries extracted from Sentinel-1 data using the algorithm developed in this study were in good agreement with in-situ measurements of lake shoreline, lake outlines delineated from the corresponding Landsat images in 2015 and lake levels for Qinghai Lake. Upon analysis, it was found that the seasonal cycles of lakes exhibited drastically different patterns across the TP. For example, large size lakes (> 100 km2) reached their peaks in August−September while lakes with areas of 50−100 km2 reached their peaks in early June−July. The peaks of seasonal cycles for endorheic lakes were more pronounced than those for exorheic lakes with flat peaks, and glacier-fed lakes with additional supplies of water exhibited delayed peaks in their seasonal cycles relative to those of non-glacier-fed lakes. Large-scale atmospheric circulation systems, such as the westerlies, Indian summer monsoon, transition in between, and East Asian summer monsoon, were also found to affect the seasonal cycles of lakes. The results of this study suggest that Sentinel-1 SAR data are a powerful tool that can be used to fill gaps in intra-annual lake observations.
ZHANG Yu, ZHANG Guoqing
The land cover dataset of Pan third pole major cities contains 14 cities (Urumqi, Xining, Lanzhou, Dhaka, Kathmandu, Lucknow, Delhi, Lahore, Islamabad, Kabul, Dushanbe, Tashkent, Bishkek and Almaty) in 2000 / 2010 / 2017, the spatial resolution of this dataset is 30 m. It includes vegetation, cultivated land, artificial surface, water body and others. Based on globeland30, mcd12q1 and globcover2009, the consistent regions were identified and retained. The inconsistent regions were reclassified by deep learning method, and the final classification results were obtained by fusing the above regions. The data has been verified by visual interpretation. The data are applied to the study of construction land dynamics and anthropogenic influence in Pan-Third Pole cities. Data type: grid. Projection mode: UTM projection.
Xin LI, LI Xin
Agricultural irrigation consumes a large amount of available freshwater resources and is the most immediate human disturbance to the natural water cycle process, with accelerated regional water cycles accompanied by cooling effects. Therefore, estimating irrigation water use (IWU) is important for exploring the impact of human activities on the natural water cycle, quantifying water resources budget, and optimizing agricultural water management. However, the current irrigation data are mainly based on the survey statistics, which is scattered and lacks uniformity, and cannot meet the demand for estimating the spatial and temporal changes of IWU. The Global Irrigation Water Use Estimation Dataset (2011-2018) is calculated by the satellite soil moisture, precipitation, vegetation index, and meteorological data (such as incoming radiation and temperature) based on the principle of soil water balance. The framework of IWU estimation in this study coupled the remotely sensed evapotranspiration process module and the data-model fusion algorithm based on differential evolution. The IWU estimates provided from this dataset have small bias at different spatial scales (e.g., regional, state/province and national) compared to traditional discrete survey statistics, such as at Chinese provinces for 2015 (bias = −3.10 km^3), at U.S. states for 2013 (bias = −0.42 km^3), and at various FAO countries (bias = −10.84 km^3). Also, the ensemble IWU estimates show lower uncertainty compared to the results derived from individual precipitation and soil moisture satellite products. The dataset is unified using a global geographic latitude and longitude grid, with associated metadata stored in corresponding NetCDF file. The spatial resolution is about 25 km, the time resolution is monthly, and the time span is 2011-2018. This dataset will help to quantitatively assess the spatial and temporal patterns of agricultural irrigation water use during the historical period and support scientific agricultural water management.
ZHANG Kun, LI Xin, ZHENG Donghai, ZHANG Ling, ZHU Gaofeng
Land surface temperature (LST) is a key parameter in the study of surface energy balance. It is widely used in the fields of meteorology, climate, hydrology, agriculture and ecology. As an important means to obtain global and regional scale LST information, satellite (thermal infrared) remote sensing is vulnerable to the influence of cloud cover and other atmospheric conditions, resulting in temporal and spatial discontinuity of LST remote sensing products, which greatly limits the application of LST remote sensing products in related research fields. The preparation of this data set is based on the empirical orthogonal function interpolation method, using Terra / Aqua MODIS surface temperature products to reconstruct the lst under ideal clear sky conditions, and then using the cumulative distribution function matching method to fuse era5 land reanalysis data to obtain the lst under all-weather conditions. This method makes full use of the spatio-temporal information of the original MODIS remote sensing products and the cloud impact information in the reanalysis data, alleviates the impact of cloud cover on LST estimation, and finally reconstructs the high-quality global 0.05 ° spatio-temporal continuous ideal clear sky and all-weather LST data set. This data set not only realizes the seamless coverage of space-time, but also has good verification accuracy. The reconstructed ideal clear sky LST data in the experimental areas of 17 land cover types in the world, the average correlation coefficient (R) is 0.971, the bias (bias) is -0.001 K to 0.049 K, and the root mean square error (RMSE) is 1.436 K to 2.688 K. The verification results of the reconstructed all-weather LST data and the measured data of ground stations: the average R is 0.895, the bias is 0.025 K to 2.599 K, and the RMSE is 4.503 K to 7.299 K. The time resolution of this data set is 4 times a day, the spatial resolution is 0.05 °, the time span is 2002-2020, and the spatial range covers the world.
ZHAO Tianjie, YU Pei
Surface soil moisture (SSM) is a crucial parameter for understanding the hydrological process of our earth surface. Passive microwave (PM) technique has long been the primary choice for estimating SSM at satellite remote sensing scales, while on the other hand, the coarse resolution (usually >~10 km) of PM observations hampers its applications at finer scales. Although quantitative studies have been proposed for downscaling satellite PM-based SSM, very few products have been available to public that meet the qualification of 1-km resolution and daily revisit cycles under all-weather conditions. In this study, therefore, we have developed one such SSM product in China with all these characteristics. The product was generated through downscaling of AMSR-E and AMSR-2 based SSM at 36-km, covering all on-orbit time of the two radiometers during 2003-2019. MODIS optical reflectance data and daily thermal infrared land surface temperature (LST) that have been gap-filled for cloudy conditions were the primary data inputs of the downscaling model, in order to achieve the “all-weather” quality for the SSM downscaling outcome. Daily images from this developed SSM product have achieved quasi-complete coverage over the country during April-September. For other months, the national coverage percentage of the developed product is also greatly improved against the original daily PM observations. We evaluated the product against in situ soil moisture measurements from over 2000 professional meteorological and soil moisture observation stations, and found the accuracy of the product is stable for all weathers from clear sky to cloudy conditions, with station averages of the unbiased RMSE ranging from 0.053 vol to 0.056 vol. Moreover, the evaluation results also show that the developed product distinctly outperforms the widely known SMAP-Sentinel (Active-Passive microwave) combined SSM product at 1-km resolution. This indicates potential important benefits that can be brought by our developed product, on improvement of futural investigations related to hydrological processes, agricultural industry, water resource and environment management.
SONG Peilin, ZHANG Yongqiang
The data set is based on a series of microwave remote sensing data, including Special Sensor Microwave Imager (SSM/I), Advanced Microwave Scanning Radiometer for Earth Observation System (AMSR-E), etc., which can be used as a reference for primary productivity. The data is from Liu et al. (2015), and the specific calculation method is shown in the article. The source data range is global, and Tibetan Plateau region is selected in this data set. This data set is often used to evaluate the temporal and spatial patterns of vegetation greenness and primary productivity, which has practical significance and theoretical value.
LIU Yi
The distribution of lakes in space and its change over time are closely related to agricultural, environmental and ecological issues, and are critical factors for human socio-economic development. In the past decades, satellite based remote sensing has been developed rapidly to provide essential data sources for monitoring temporal lakes dynamics with its advantage of rapidness, wide coverage, and lower cost. This dataset was produced from Landsat images using the automated water detection method (Feng et al, 2015). We collected 96,278 Landsat images (about 25 terabytes) that acquired since 2000 with less than 80% cloud contamination in the arid region of central Asia and Tibetan Plateau. Water is detected in each of the image and then aggregated to monthly temporal resolution by taking advantage of the high-performance processing capability and large data storage provided by Global Land Cover Facility (GLCF) at University of Maryland. The results are validated systematically and quantitatively using manually interpreted dataset, which consists of a set of locations collected by a stratified random sampling strategy to effectively represent different spatial-temporal distributions in the region. The validation suggests high accuracy of the results (overall accuracy: 99.45(±0.59); user accuracy: 85.37%±(3.74); produce accuracy: 98.17(±1.05)).
FENG Min, CHE Xianghong
The dataset include ground-based passive microwave brightness temperature, multi-angle brightness temperature, ten-minute 4-component radiation and snow temperature, daily snow pit data and hourly meteorological data observed at Altay base station（lon：88.07、lat: 44.73）from November 27, 2015 to March 26, 2016. Daily snow pit parameters include: snow stratification, stratification thickness, density, particle size, temperature. These data are stored in five NetCDF files: TBdata. nc, TBdata-multiangle. nc, ten-minute 4 component radiation and snow temperature. nc, hourly meteorological and soil data. nc and daily snow pit data.nc. TBdata. nc is brightness temperature at 3 channels for both polarizations automatically collected by a six-channel dual polarized microwave radiometer RPG-6CH-DP. The contents include Year, month, day, hour, minute, second, Tb1h, Tb1v, Tb18h, Tb18v, Tb36h, Tb36v, incidence angle, azimuth angle. TBdata-multiangle.nc is 7 groups of multi-angle brightness temperatures at 3 channels for both polarizations. The contents include Year, month, day, hour, minute, second, Tb1h, Tb1v, Tb18h, Tb18v, Tb36h, Tb36v, incidence angle, azimuth angle. The ten-minute 4 component radiation and snow temperature.nc contains 4 component radiation and layered snow temperatures. The contents include Year, month, day, hour, minute, SR_DOWN, SR_UP, LR_DOWN, LR_UP, T_Sensor, ST_0cm, ST_5cm, ST_15cm, ST_25cm, ST_35cm, ST_45cm, ST_55cm. The hourly meteorological and soil data.nc contains hourly weather data and layered soil data. The contents include Year, month, day, hour, Tair, Wair, Pair, Win, SM_10cm, SM_20cm, Tsoil_5cm, Tsoil_10cm, Tsoil_15 cm, Tsoil_20cm. The daily snow pit data.nc. is manual snow pit data. The observation time was 8:00-10:100 am local time. The contents include Year, month, day, snow depth, thickness_layer1, thickness_layer2, thickness_layer3, thickness_layer4, thickness_layer5, thickness_layer6, Long_layer1, Short_layer1, Long_layer2, Short_layer2, Long_layer3, Short_layer3, Long_layer 4, Short_layer4, Long_layer5, Short_layer5, Long_layer6, Short_layer 6, Stube, Snow shovel_0-10, Snow shovel _10-20, Snow shovel _20-30, Snow shovel _30-40, Snow shovel _40-50, Snow fork_5, Snow fork _10, Snow fork _15, Snow fork_20, Snow fork_25, Snow fork_30, Snow fork_35, Snow fork_40, Snow fork_45, Snow fork_50, shape1, shape2, shape3, shape4, shape5,
DAI Liyun
Based on AVHRR-CDR SR products, a daily cloud-free snow cover extent dataset with a spatial resolution of 5 km from 1981 to 2019 was prepared by using decision tree classification method. Each HDF4 file contains 18 data elements, including data value, data start date, longitude and latitude, etc. At the same time, to quickly preview the snow distribution, the daily file contains the snow area thumbnail, which is stored in JPG format. This data set will be continuously supplemented and improved according to the real-time satellite remote sensing data and algorithm update (up to may 2019), and will be fully open and shared.
HAO Xiaohua
This dataset includes one scene acquired on (yy-mm-dd) 2012-07-25, covering the natural oasis eco-hydrology experimental area in the lower reaches of the Heihe River Basin. This datum contains panchromatic and multi-spectral bands, with spatial resolution of 0.6 m and 2.4 m, respectively. The data product level of this image is Level 2A. QuickBird dataset was acquired through purchase.
LI Xin
On 19 August 2012, a Leica ALS70 airborne laser scanner boarded by the Y-12 aircraft was used to obtain the point cloud data. Leica ALS70 airborne laser scanner has unlimited numbers of returns intensities measurements including the first, second, third return intensities. The wavelength of laser light is 1064 nm. The absolute flight altitude is 2900 m with the point cloud density 1 point per square meter. Aerial LiDAR-DSM was obtained through parameter calibration, automatic classification of point cloud density and manual editing.
XIAO Qing, Wen Jianguang
This is the vegetation index (NDVI) for Maduo County in July, August and September of 2016. It is obtained through calculation based on the multispectral data of GF-1. The spatial resolution is 16 m. The GF-1 data are processed by mosaicking, projection coordinating, data subsetting and other methods. The maximum synthesis is then conducted every month in July, August, and September.
LI Fei, Fei Li, Zhijun Zhang
Data content: precipitation data of the Aral Sea basin from 2015 to 2018. Data sources and processing methods: from the new generation of global precipitation measurement (GPM) of NASA (version 06, global precipitation observation program), the daily rainfall can be obtained by adding the three-hour rainfall data, and then the eight day rainfall can be obtained. Data quality: the spatial resolution is 0.1 ° x 0.1 ° and the temporal resolution is 8 days. The value of each pixel is the sum of rainfall in 8 days. Data application results: under the background of climate change, it can be used to analyze the correlation between meteorological elements and vegetation characteristics.
XIAO Qing, Wen Jianguang
Based on a large number of measured aboveground biomass data of grassland, the temperate grassland types were divided according to the vegetation type map of China in 1980s Based on the Landsat remote sensing data of engine platform, the random forest model of grassland aboveground biomass and remote sensing data was constructed for different grassland types. On the basis of reliable verification, the annual aboveground biomass of grassland from 1993 to 2019 was estimated, and the annual spatial data set of aboveground biomass of temperate grassland in Northern China from 1993 to 2019 was formed. Aboveground biomass is defined as the total amount of organic matter of vegetation living above the ground in unit area. The original grid value has been multiplied by a factor of 100, unit: 0.01 g / m2 (g / m2). This data set can provide a scientific basis for the dynamic monitoring and evaluation of temperate grassland resources and ecological environment in northern China.
ZHANG Na
The global Cryosat-2 GDR dataset is generated by the European Space Agency (ESA); it has a temporal coverage from 2010 to 2016 and covers the globe. On April 8, 2010, the ESA launched the Cryosat-2 high-tilt polar orbit satellite. The satellite is equipped with an SAR Interferometer Radar Altimeter (SIRAL), which is mainly used to monitor polar ice thickness and sea ice thickness changes, and, furthermore, to study the effects of melting polar ice on global sea level rise and that of global climate change on Antarctic ice thickness. The altimeter operates in the Ku-band and at a frequency of 13.575 GHz, it includes three measurement modes. One is a low-resolution altimeter measurement mode (LRM) that points to the subsatellite point to obtain all surface observations for land, sea, and ice sheets; its processing is similar to ENVISAT/RA-2, with an orbital resolution of 5 to 7 km. The second is the Synthetic Aperture Radar (SAR) measurement mode, which is mainly used to improve the accuracy and resolution of sea ice observations; it can make the resolution along the orbit reach approximately 250 m. The third is the Interferometric Synthetic Aperture Radar (InSAR), which is mainly used to improve the accuracy of areas with complex terrain such as the edges of ice sheets or ice shelves. The CryoSat -2/SIRAL data products mainly include 0-level data, 1b-level data, 2-level data and high-level data. The Cryosat-2/SIRAL products consist of two files: an XML head file (.HDR) and a data product file (.DBL). The HDR file is an auxiliary ASCII file for fast identification and retrieval of the data files. 1b-level products are stored separately according to the measurement modes, and the data recording formats of different modes are also different. Each waveform in LRM mode and SAR mode has 128 sampling points, while that in SARIn mode has 512 sampling points. 2-level GDR products are available for most scientific applications, including measurement time, geographic location, altitude, and more. In addition, the altitude information in GDR products has been obtained through instrumental calibration, transmission delay corrections, geometric corrections, and geophysical corrections (such as atmospheric corrections and tidal corrections). The GDR products are single global full-track data, that is, the measurement results of the three modes. After different processing, they are combined in chronological order; thereby, the data recording formats are unified. The data in the three modes use different waveform retracking algorithms to obtain altitude values. In the latest updated Baseline C data, the LRM mode data use three algorithms: Refined CFI, UCL and Refined OCOG.
SHEN Guozhuang, FU Wenxue
This dataset includes component temperatures measured by the thermal infrared (TIR) radiometers at the Mixed Forest and Sidaoqiao stations between 22 July, 2014 and 19 July, 2016. The Mixed Forest (101.1335 °E, 41.9903 °N, 874 m.a.s.l.) and Sidaoqiao (101.1374 °E, 42.0012 °N, 873 m.a.s.l.) stations were located in the downstream of the Heihe River basin, Dalaihubu Town, Ejin Banner, Inner Mongolia. At the Mixed Forest station, two TIR radiometers (SI-111, Apogee Instruments Inc., USA) connected to a data logger (CR800, Campbell Scientific Inc., USA) measured component temperatures of the sunlit canopy and shaded canopy. TIR radiometers were mounted horizontally at 5 m height on iron rods just south and north of a tree and pointed to its canopy. The distance from the sensor to the canopy was ~1 m. At the Sidaoqiao station, two SI-111 TIR radiometers connected to a CR800 data logger measured component temperatures of the soil and shrub. The first sensor pointed from 2 m height under a viewing zenith angle of 45° to bare soil; the second sensor was mounted at 1-m height and pointed horizontally into the shrub canopy.
ZHOU Ji, LI Mingsong , MA Jin
This dataset includes component temperatures measured by the thermal imager at the Mixed Forest and Sidaoqiao stations between 23 July and 18 August, 2014. The Mixed Forest (101.1335 °E, 41.9903 °N, 874 m.a.s.l.) and Sidaoqiao (101.1374 °E, 42.0012 °N, 873 m.a.s.l.) stations were located in the downstream of the Heihe River basin, Dalaihubu Town, Ejin Banner, Inner Mongolia. At the Mixed Forest station, a Testo 890-2 thermal imager (Testo Inc., Germany) with a resolution of 640 × 480 pixels was employed to acquire brightness temperature images. The imager was manually operated from a 10-m height platform of the tower between 10:00-16:00 (China Standard Time, CST) with an observation interval of 1-h on cloudless days. On August 4th observations were acquired between 11:00 and 17:00 at an interval of 10-min to match observations acquired with an airborne TIR imager. The ground based imager was pointed to five viewing directions (southeast-SE, east-E, northeast-NE, northwest-NW, and southwest-SW) and was inclined 25°–45° below the horizon depending on viewing direction. At Sidaoqiao station, a Testo 875-2i imager (Testo Inc., Germany) with a resolution of 160 × 120 pixels was manually operated from a 10-m high platform to acquire brightness temperature images in directions SW, SE, NE, and NW. Depending on the targets in each viewing direction, the imager was inclined to 30°–45° below the horizon. Observations at Sidaoqiao and Mixed Forest stations were almost synchronous. Furthermore, visible images were taken simultaneously with the aforementioned two TIR imagers (2048 × 1536 pixels for Testo 890-2 and 640 × 480 pixels for Testo 875-2i).
LI Mingsong , MA Jin
This dataset includes three scenes, covering the artificial oasis eco-hydrology experimental area of the Heihe River Basin, which were acquired on (yy-mm-dd hh:mm, BJT) 2012-07-25 07:12, 2012-07-28 19:55, 2012-08-02 07:12. The data were all acquired at PingPong mode with product level of SLC, and these three images are of VV/VH, HH/HV and VV/VH polarization, respectively. COSMO-SkyMed dataset was acquired from Italian Space Agency (ASI) “COSMO-SkyMed project 1720: HYDROCOSMO” (Courtesy: Prof. Shi Jiancheng from the State Key Laboratory of Remote Sensing Science of China).
Agenzia Spaziale Italiana (ASI)
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
